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Abstract. - The disordering of an initially phase segregated system of finite size, induced by 
the presence of highly mobile vacancies, is shown to exhibit dynamic scaling in its late stages. A 
set of characteristic exponents is introduced and computed analytically, in excellent agreement 
with Monte Carlo data. In particular, the characteristic time scale, controlling the crossover 
between increasing disorder and saturation, is found to depend on the exponent scaling the 
number of vacancies in the sample. 



The structure and formation of surfaces and interfaces is of great scientific and technical 
importance and has attracted considerable interest over the past decade [Q. A second problem 
of similar relevance concerns the morphology and dynamics of ordering in bulk systems after 
a rapid temperature quench|^. The evolution of spatial structures in both of these processes 
exhibits dynamic scaling in the late-time regime so that typical configurations at different times 
are self-similar after an appropriate rescaling of space and time. Clearly, the "inverse" scenario, 
i.e., the destruction of interfaces and the bulk disordering of an initially phase-segregated 
system, are also of major significance, being related to natural erosion phenomena . Mobile 
defects, present in many materials, can obviously play a major role in the disordering process, 
especially if their dynamics is fast on the time scale governing the bulk particles. In this Letter, 
we consider a simple model for defect-mediated interface destruction and bulk disordering in a 
finite system, corresponding to a real material in which the characteristic time scale for vacancy 
diffusion is much faster than the ordinary bulk diffusion time. A fixed number of vacancies are 
initially located at a smooth, nonelastic interface separating two ordered bulk phases. These 
vacancies perform a homogeneous Brownian random walk, by exchanging positions with the 
otherwise completely passive bulk particles. As a result, the interface is gradually destroyed, 
and the system approaches a completely disordered equilibrium state. Three distinct temporal 
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regimes, separated by two crossover times, are observed. Our key result is that the late stages of 
this process exhibit dynamic scaling. The characteristic exponents are computed analytically. 
In particular, the late crossover time is found to scale as a power of system size, determined 
by the spatial distribution of the defects. 

This Letter is organized as follows: We first describe our model, followed by a summary of 
our simulation results. Turning to analytic methods, we characterize the initial configuration 
and the final state exactly. The intermediate regime is described by a set of mean-field 
equations of motion which correctly predict both the scaling form for an appropriately defined 
"disorder parameter" and the characteristic exponents controlling its size and time dependence. 
We conclude with some comments and questions. More details will be published elsewhere 
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Fig. 1. - An initially smooth interface (a) gradually becomes rough (b) (2 • 10'^ steps), and 
destroyed after (c) (after 5 • 10^ steps), (d) represents the completely mixed, equilibrium state 
(more than 2 • 10^ steps). Here, L — 60. 



Our model is defined on a two-dimensional square lattice of sites, with L ranging from 10 
to 100. Each lattice site can be occupied either by a white particle, a black particle or a hole 
(vacancy), so that multiple occupancy is forbidden. The number of vacancies is denoted by M . 
For simplicity, we restrict our attention only to the case of L^/2 black particles and L^/2 — M 
white particles. Since we wish to describe bulk disordering rather than the "evaporation" of the 
two bulk phases into an empty region, we will keep the number of holes small, i.e., M/L^ << 1. 
Initially, the upper (lower) half of the system is occupied by white (black) particles only, and 
the vacancies are randomly distributed along the flat interface. At each Monte Carlo step 
(MCS), a hole is chosen at random and changes place, with equal probability, with one of its 
four nearest neighbors. All particle-particle exchanges are forbidden. The boundary conditions 
at the right and left edge are periodic; for the top and bottom edges, we have simulated both 
brickwall and periodic boundary conditions, without detecting any significant differences in 
the scaling form and exponents. 

We are interested in the dynamics of the interface destruction and bulk disordering caused 
by these vacancies. Monitoring the evolution of a typical configuration in a 60x60 lattice, shown 
in Fig.l, we observe that the interface is gradually destroyed, via increasing fragmentation, 
until it ceases to be identifiable. As a quantitative measure for the growing disorder, we 
consider the total "surface area" of the black regions, i.e., the total num&er of black- white and 
black- hole bonds, after t MCS. In analogy with the Ising model, we will refer to these pairs 
as "broken bonds" . Averaging this quantity over many configurations (200 for L < AO and 20 
otherwise) we obtain the "disorder parameter", A{L,AI,t), which depends on system size L, 
vacancy number M, and 'time' t. Since M / ^ 1, the dominant contribution to this quantity 
is just the average number of black-white bonds. Fig. 2a shows yl(40, l,t), i.e., the average 
surface area in a 40 x 40 system containing a single vacancy. One clearly distinguishes three 
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regimes, drawn schematically in the inset: an early regime (I), the intermediate, or scaling, 
regime (II), and finally a late or saturation regime (III) in which the system has effectively 
reached the steady state. Physically, the three regimes are easily interpreted: Tracking the 
position of the vacancy, we observe that, for early times, it is confined to a small (yet growing) 
region centered on its starting point, far from the boundaries of the system. Thus, the early 
regime is strongly dependent on the dimensionality of the system (see |q|). 
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Fig. 2. - (a) the total number of broken bonds, A{L, M, t) vs time t (i.e., the number of steps) 
for the single vacancy case M = 1 and system size L = 40. (b) scaling plot of A{L, 1, t)/L'^ vs 
t/L'^, for system sizes ranging from i = 10 to L = 40. 



In contrast, for intermediate and late times, the vacancy is found with essentially equal 
probability anywhere in the system. The difference between intermediate and late stages lies 
in the degree of disorder among the black and white particles. 

Our key observation is that, independent of dimension. Regimes II and III exhibit dynamic 
scaling. To characterize this behavior, we define a set of exponents: First, the saturation value 
of A scales with system size according to limf^oo A{L, 1, t) = Asat{L, 1) ~ L". Second, in the 
intermediate regime, A{L, 1, t) grows as A{L, l,t) ^ L°'t^ . The L-dependent factor is required 
here to absorb a nontrivial shift along the hit axis in dimensions other than two in order to 
collapse curves associated with different system sizes. Finally, the two crossover times ("early" 
and "late") scale as te ~ and ti ~ L^'. In d = 2, the data give a = 2 ±0.1, f3 = 0.5 ±0.06, 
Ze — 2zL 0.2, z; = 4 ± 0.2, a = ± 0.1. Since a is effectively zero, excellent data collapse is 
obtained by plotting A{L, 1, t)/L°' versus t/L^\ shown in Fig. 2b for a range of system sizes. 
Thus, the intermediate-to-late time crossover can be summarized by the usual Family- Vicsek 
[H scaling form 

A{L,l,t) ^ (1) 

with a scaling function /, satisfying f{x) ~ const for a: ^ 1, and f{x) ^ for a; <C 1 
(but large enough to fall within Regime II). The consistency condition A{L, 1, ti) ~ Asat{L, 1) 
yields the scaling law 



cr ± = a, 



(2) 
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which is manifestly satisfied by the measured exponents. The same indices are observed for 
several vacancies, provided their number remains constant as the system size L is varied. We 
stress that, due to the presence of the novel exponent a, our scaling law (|^) is distinct from 
the familiar zP = a which controls surface growth in, e.g., the Edwards- Wilkinson or KPZ 
[| models. 

In a real system it is to be expected that the number of vacancies itself depends on the 
system size. For generality, we allow M ^ U' , where 7 G [0, d] will be called the vacancy 
number exponent. The data reported above then correspond to 7 = (and d — 2). Another 
natural case is 7 = d — 1, which could arrise if we imagine welding two different solids along 
a common surface with vacancies uniformly distributed thereupon. If the two surfaces were 
of equal fractal dimensions, a uniform distribution of vacancies on the interface would lead to 
a non- integer 7. In this study wc have examined the case 7 = 1, i.e., the vacancy number 
being extensive in the initial interfacial length, M oc L. Here, the data collapse only if (i) 
curves corresponding to different system sizes are shifted by an M-dependent factor along the 
Ini-axis, and (ii) the late crossover exponent takes on the new value = 3 ± 0.2, while a 
and /3 remain unchanged! This behavior suggests that the exponent a acquires an additional 
significance, depending non-trivially on 7. For 7 = 1, the data are consistent with a = 1/2, 
so that the scaling form ([|) still holds and the scaling law (||) is satisfied with the modified 
exponent z;. 

In order to understand the role of 7, we now turn to an analytic description of the disordering 
process, in general dimension d. For simplicity, we choose fully periodic boundary conditions 
on a hypercube of side length L. The initial value, A{L, M, i = 0) = 2L'^~^, reflects two fiat 
interfaces between black and white particles, with vacancies replacing M of the white particles 
at one of the interfaces. Similarly, we have Asat{L, M) = |L^''/(L'' — 1) in the saturation 
regime, since the final steady state is completely random. Thus, we can read ofi' the exact 
result a = d. 

The remainder of our analysis will be a mean-field theory, based on a set of equations of 
motion for the coarse-grained local hole and particle densities [|[ Q which are most conveniently 
discussed in continuous time and space, in general dimension d. A well-defined continuum limit 
is obtained by letting the lattice constant a vanish and identifying the microscopic time scale 
with T = jld. In this limit, physical distances such as the system size L remain fixed while 
the number of sites approaches infinity. For convenience, we define V = L'^ . 

First, we note that each of the vacancies performs a simple random walk. As a result, the 
evolution equation for 0(x, <), the density of vacancies, at position x = (xi, Xd) and time t is 
just a diffusion equation, 9t0(x, i) = V^(/i(x, t) , subject to fully periodic boundary conditions 
on (/)(x, t) and some initial condition. Choosing the normalization <?!>(x, t) = M , we see that, 
for our simulations, (f>{x, 0) — j^L 1 5 {y — L/2) where y = Xd- Note that we still use the symbol 
t here, since the Monte Carlo time and its continuum limit differ only by the scale factor r. 
The final state is, of course, trivial: 0(x, 00) = M/V. 

A similar evolution equation for the black particle density, -(/'(x, t), is easily derived from a 
microscopic master equation. Truncating all correlations, we find 

a-'*5t^(x, t) = 0(x, i) V'V^x, t) - ^/.(x, t)VV(x, t). (3) 

This equation simply tallies the local change in black particle density at position x: the first 
term reflects a gain, provided a vacancy is initially present and a black particle "diffuses" 
in from a neighboring site. The second term accounts for a loss, due to a black particle 
jumping from x to a vacant nearest-neighbor site. Similar equations have been discussed in 
the context of biased diffusion of two species g. The prefactor a appears for dimensional 
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reasons. Choosing Jy t) — V/"^, we write the initial condition for the half-filled system as 
V'(x, 0) = 1 — 9{y — L/2). Of course, the final state is given by oo) = 1/2. 

The key difference between these two equations resides in a separation of time scales: the 
hole density relaxes very fast, driven by a diffusion coefficient of unity, compared to the particle 
density whose diffusion "coefficient" is (/)(x, t) = 0{1/V). Performing a simple random walk, 
the vacancies reach the edge of the system after a time of order L^. Thus, we identify the 
early crossover time te oc and read off Ze — 2. For later times, the holes are uniformly 
distributed over the system, and we may replace by its final value, so that reduces 
to a simple diffusion equation, dtip = DW^ip, where D = a'^M/V. Its solution, subject to the 
initial and fully periodic boundary conditions, is easily found: 



7/;(x,<) 



E 



sin[2n{2n-l)y/L] ,t{2n~if 



2n-l 



(4) 



where e = An^D/L'^. 

To make contact with the "disorder parameter" , A, measured in the simulations, we could 
define an operator on the lattice configurations for the total number of broken bonds, make 
a meanfield approximation for pair correlations, and take the continuum limit jij. Here, 
let us present a short cut, based on the analogy with the Ising model and its associated 
coarsed grained version. The local magnetization density, $(x), of the former clearly maps 
into ip{x.,t) — i for our case. Meanwhile, the local energy density of the Ising model is given 
precisely by the broken bonds. In the standard literature |I0|, the total energy of the Ising 
model is written as Ki — K2 Jy where the constants Ki may be fixed by the scale and the 
"zero" of the energy. Since our interest is also the total number of broken bonds, we may write 
A oc KV ~ Jy [■(/'(x, ^) ~ 5] • Note that this form expresses the extensivity of A. Of course, 
since the total number of broken bonds at i = is only ©(L'^"^), we set the initial value of A 
to be 0. Demanding the final value be dV/2, we have 



A{L,M,t) 



-V ~2d 



V'(x,t) 



Using (^, we find the time evolution of this quantity: 



(5) 



AiL,M,t) = -V[l-S{2et)] 



(6) 



where 5(0 = ^ e"^*^""^^V(2f^ - 1)^- Since S{0) = I and S{oo) = 0, we verify that 
A does take on the correct initial and final values. To capture the time dependence in the 
intermediate regime (II), i.e., et <C I, we reexpress the infinite sum via a Poisson resummation 
[ pr] . Introducing Um = 7rm/2-\/2ei, we find 



A{L,M, t) 



2d 



VV2et{ 1 + 



00 

E 



(-ir 



(7) 



where r(»,») denotes the incomplete Gamma function. In this form, the sum over m is 
suppressed for small et. Thus, A{L,M,t) cx VV2ei oc L'^^/MtJT^, yielding the remaining 
exponents, namely (3 — ^ independent of dimension, and a — ^{d + j — 2). The late crossover 
time, naturally defined by eti = I, scales as ti ^ ip'+d-^i whence one obtains zi = 2 + d — j. 
In two dimensions our theory predicts = 4 and a = for 7 = 0, while zi — 3 and a 
7 = 1 , in complete agreement with the Monte Carlo data. 



ifor 



6 



EUROPHYSICS LETTERS 



To summarize, we have analyzed the disordering process of an initiaUy phase segregated 
system, driven by highly mobile Brownian vacancies distributed according to the exponent 7. 
The late stages of the evolution exhibit dynamic scaling. A set of exponents {ze, z;, a, cr; /3} 
can be defined, characterizing, respectively, the system size dependence of two crossover times, 
the final saturation value of the number of broken bonds, A{L, M, t), and its amplitude in the 
intermediate regime. The temporal growth of A during the latter regime is captured by the 
exponent /?. All indices can be computed analytically, in excellent agreement with the data. 
Our key result is that the typical time scale ti ~ L^', which controls the crossover between 
increasing disorder and saturation, is set by z; — 2 + d — 7, and thus depends explicitly on 
both the space and fractal dimensionalities, d and 7. Measurements of z; can therefore provide 
information about the vacancy distribution in a sample. In the most familiar case, standard 
vacancy diffusion in solids, the number of vacancies is extensive (7 = d), so that the well-known 
result = 2 is reproduced Even though our model is very simple, it forms the basis for 
the description of a large variety of related problems. Work is in progress to analyze the effect 
of external driving forces, interparticle interactions and vacancy-induced catalytic reactions. 
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